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ABSTRACT 

We study the effect of surface brightness on the mass-metallicity relation using nearby galaxies whose gas 
content and metallicity prohles are available. Previous studies using fiber spectra indicated that lower surface 
brightness galaxies have systematically lower metallicity for their stellar mass, but the results were uncertain 
because of aperture effect. With stellar masses and surface brightnesses measured at WISE W1 and W2 bands, 
we re-investigate the surface brightness dependence with spatially-resolved metallicity prohles and hnd the 
similar result. We further demonstrate that the systematical difference cannot be explained by the gas content 
of galaxies. For two galaxies with similar stellar and gas masses, the one with lower surface brightness tends to 
have lower metallicity. Using chemical evolution models, we investigate the inhow and outhow properties of 
galaxies of different masses and surface brightnesses. We hnd that, on average, high mass galaxies have lower 
inhow and outhow rates relative to star formation rate. On the other hand, lower surface brightness galaxies 
experience stronger inhow than higher surface brightness galaxies of similar mass. The surface brightness 
effect is more signihcant for low mass galaxies. We discuss implications on the different inhow properties 
between low and high surface brightness galaxies, including star formation efficiency, environment and mass 
assembly history. 

Subject headings: galaxies: abundances, galaxies: spiral, galaxies: formation 


1. INTRODUCTION 

Chemical enrichment is one of the keys to understand the 
evolution of galaxies. Heavy elements, or metals, are synthe¬ 
sized in stars then released into the interstellar medium (ISM) 
by stellar winds or supernova explosions. The metal-enriched 
ISM subsequently acts as the raw material for the next gen¬ 
eration of new-born stars. Tight relations between the stellar 
mass and both gas-phase and stellar metallicities were discov¬ 
ered, where galaxie s with larger stellar mass on average have 
higher metallicities ^remonti et al.|2004[|Gallazzi et al.|2005[ 


Savaglio et al.|2005||Lee et al. 12006)1. 

To hrst order, the origin of the mass-metallicity relation can 
be understood as a simple process of recycling of metals in the 
ISM. Galaxies with larger stellar mass have synthesized and 
released more heavy element throughout their lives, therefore 
higher metallicity is naturally expected. Nevertheless, this 
simple picture does no t tell the whole story. As pointed out by 
Tremonti et al.|(|2004|), the dispersion in the mass-metallicity 
relation cannot be explained simply by measurement uncer¬ 
tainties of metallicity. There must be other factors affecting 
the metallicity. 

Several studies have been searching for parameters beyond 
stellar mass which affects metallicity. An anti-correlation be¬ 
tween gas-phase metallicity and star formation rate (SFR) has 
been found both locally and at high redshifts, where at hxed 
stellar mass, galaxies with higher SFR have lower metallicity 
( [Ellison et al.|2008[JMannucci et al.|2010[ [Lara-Lopez et al.j 

Z010||Yabe et al.|2015|l. A similar anti-correlation seems also 


to be present between the gas content and metallic ity, where 
gas-rich galaxies have on average lower metallicity (jBothwell 
et al.|20l3| Hughes et al.|2013l. There is also a dependency 
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on galaxy structure, where larger, or lower surface brightness, 
galaxies tend to have lower metallicity for their stellar mass 


( Effison et al.|200^ Liang et al.|2010| Salim et al.|2014 1 . 
Most of these results are obtained with data from SI 


SDSS 

hber spectra. Although the sample size is large and pro¬ 
vides good statistics, data from hber spectra have fundamen¬ 
tal issues. The spectroscopic information, including SFR and 
metallicity, comes from only the center of the galaxy, while 
stellar mass, gas mass and galaxy structure are derived from 
imaging with different apertures. Moreover, the coverage 
varies among galaxies, depending on the size, structure and 
distance of galaxies. Therefore, some of the dependencies on 


a second parameter could come from aperture effects ( Ellison 
et al.||20()8| [Hughes et al.||2013[ [Sanchez et al.[[2013 1 . This 
calls for an approach, where spatially resolved information of 
metallicity is used. 

In the interpretation of the observed relationship between 
integrated metallicity and stellar mass, it has become clear 
that a simple closed box chemical evolution model with no 
exchange of material between the galaxies and their outside 
environment ([Searle & Sargent|1972[[Pagel & Patchett|1975|l 
is insufficient. In reality, galaxies experience gas inhows from 
mergers and accretion and outhows caused by feedback from 
supernova, starburst or AGN and, thus the close-box assump¬ 
tion does not hold. Exchanging material between galaxies and 
the environment modihes the metallicity of galaxies. By in¬ 
troducing inflow and outflow activities, theoretical models are 
able to reproduce qu antitatively the observed integrated mass- 
metallicity relation (|Spitoni et al.||2010[ Peeples & Shank^ 


[20^[Zahid et al.|2012[|Lilly et al.|2013[|Zahid et al.|2014| l. 

With spatially resolved observations, to overcome the un- 


certainties mentioned above a modified theoretical approach 
is also required, which utilizes the additional information 
coming from measurement of radially resolved metallicity 
profiles . Recently, a variety of such models has been pre¬ 
sented (Pilkington et al.|2012[[Mott et al.|20T3j[Kudritzki et[ 


sentea (uiiKingion et ai.|zuizniviott et ai.|zui5HJS.uariizKi et| 
[al.[2015|l. This allows to impose constraints on the inflow and 
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outflow pro perties of individual galaxies by chemical e volu- 
tion models ( jAscasibar et al.|2015[|Kudritzki et al.|2015| l, and 
then to investigate whether there is any dependence on other 
physical properties of galaxies. 

In this paper, we focus on the surface brightness effect on 
the mass-metallicity relation. The surface brightness, or sur¬ 
face density, of a galaxy de pends on the angular momentum 
of galaxy for a given mass ( Dalcanton et al.|ri997||Mo et al. 
1998| l, whic h is affected by the mass loss an d accretion history 
of a galaxy ( |Dutton & van den Bosch|2012| l. We will examine 
the interplay among surface brightness, metallicity, and prop¬ 
erties of inflow and outflow of galaxies from nearby galaxies, 
where the spatially-resolved information is available. In Sec¬ 
tion 2, we describe the data sets used in the analysis. Section 3 
presents the dependence of metallicity on surface brightness 
inferred from both spectroscopy and broadband color in two 
independent samples. We use chemical evolution models to 
constrain the inflow and outflow of galaxies in Section 4, and 
discuss possible implications and systematics in Section 5. 
The summary is given in Section 6. We use AB magnitudes 
in this paper unless noted. 

2. DATA AND ANALYSIS 
2.1. Samples 

In this paper, we present two samples of galaxies with 
metallicities estimated from spectra of HII regions and broad¬ 
band colors respectively. We use data fro m surveys including 
the Wide-Field Infrared Explore (WIS E; [Wright et al.|2010| l, 
the Sloan Digital Sky Survey (SDSS; York et aTI 200011, "and 
the Arecibo Legacy East ALEA Survey (ALEALEA; [Haynes 
et al.pOlTj ). The synergy among these large surveys allows 
us study large number of galaxies with homogeneous obser¬ 
vations. The distribution of the stellar mass, gas mass, gas 
fraction, and surface brightness of the sample is presented in 
Figure [1] 

2.1.1. Spectroscopic Sample 

The first sample consists of 118 galaxies with gas-phase 
metallic ities measured from El II regions in [Pilyugin et al.[ 
([2014a[l. This paper compiled published spectra of El II re- 
gions in 130 nearby galaxies and determined the gas-phase 
metalli city of each H II region in a uniform way. For each 
galaxy, Pilyugin et al. (2014ai reported the central metallic¬ 
ity and metallicity gradient as a function of radius. This set 
of galaxies provides a homogenized sample with spatially re¬ 
solved information. We select galaxies whose HI flux is avail¬ 
able for further discussion and refer this sample as the spec¬ 
troscopic sample 

The first row of Figure [T] shows the properties of the spec¬ 
troscopic sample. The majority of galaxies in this sample have 
stellar mass M* > IO'^Mq and relatively high surface bright¬ 
ness, while the distributions have tails towards lower mass 
and lower surface brightness end. This sample covers a spe¬ 
cific range of stellar-to-gas mass ratio, 0 < log(M*/Mg) < 1. 
Few galaxies are located beyond this range. 


2.1.2. Photometric Sample 

The s econd sample is dra wn from the ALFALFA a.40 data 
release (IHaynes et al.|201 l|l. Our sample selection starts from 
ALFALFA sources with heliocentric velocity Vhei < 3000 km 
s“'. The velocity, or equivalent, distance limit is chosen as 
a compromise between sample size and completeness. At 
larger distances, only Hl-massive galaxies will be detected by 


ALFALFA. Also, the angular resolution of the WISE images 
prevents us from measuring the structure parameters of stellar 
disks for fainter or less extended galaxies, thus causing poten¬ 
tial bias. 

We then measure the magnitude and surface brightness pro¬ 
file of each ALEALEA source within the velocity limit in the 
WISE W1 3 .4rim band. The procedure will be described in 
Section 2.2 We select galaxies with m^i < 16, b > 0.8', 
where b is the minor axis that measured at = 25 mag 
arcsec^, and axis ratio b/a > 0.35. We consider measure¬ 
ments from galaxies fainter or smaller than these limits are 
not reliable. The limit on axis ratio removes edge-on galax¬ 
ies, where correction of the geometric effect of the inclination 
is less certain. 

Lastly, we excluded galaxies not covered by SDSS imaging. 
We end up with a sample of 488 galaxies. This sample cov¬ 
ers wider ranges of stellar mass, gas mass, stellar-to-gas mass 
ratio, and surface brightness than the spectroscopic sample, 
especially for the low mass, low surface brightness and gas- 
rich end of the distribution (Ligure[2. We refer to this sample 
as the photometric sample. 

2.2. Photometry and Surface Brightness Profile 

Lor this study we perform surface photometry on the WISE 
W1 and W2 band, and SDSS g', and z' band images of each 
galaxy, using elliptical apertures with fixed centers, orienta¬ 
tions and shapes with varying major axes in step of 3". Lore- 
ground stars and contaminating neighbor galaxies are masked 
and excluded from measurements. 

We determine ellipses for surface photometry from W1 im¬ 
ages as following. We use the shape from the RC3, or SDSS 
r'-band when needed, to measure the sky level for the first 
pass. We then fit an ellipse to the isophot at 25 mag arcsec“^ of 
the sky-subtracted image. If the shape and the orientation of 
the output ellipses are consistent with those of input ellipses, 
we use it as our input for surface photometry. If the shapes of 
input and output ellipses are inconsistent (AP.A. > 20 deg or 
Ah/a > 0.2), we choose the better one based on visual inspec¬ 
tion. The chosen orientation and shape are then applied for 
surface photometry at all bands. We correct the Galactic ex- 
tinction base on the extinction map ofjSchlafly & Linkbeiner 
( 201 1[ ) and estimate the internal extinction following the pre- 
scription of[Verheijen & Sancisi[([2001|l. Throughout this pa¬ 
per, we convert extinction am ong different filte rs using the 
conversion factors provided byjYuan et al.[([2013|l. 

We use the central surface brightness of the disk in the W1 
band as the surface brightness discussed in this study. The 
3.4/im W1 band is nearly free from concern for internal ex¬ 
tinction. We fit the annuli-averaged surface brightness profile 
in W1 using an exponential form of 


ji(r) = /To X exp(-r//i), 


( 1 ) 


, where /io is the central surface brightness of the disk and h is 
the scale length. To exclude the effect of the bulge, we adopt 
a fiducial fitting range between effective radius, r^., and the 
radius of isophot of 25.5 mag arcsec“^, r 25 , 5 . In some cases, 
the surface brightness at is still affected by the bulge or 
bar. We then manually adjust the inner fitting range to avoid 
the effect from s tructures in the inner pa rt of the galaxy. As 
demonstrated by [McDonald et al.[ ([2009[), we expect little to 
no statistical difference between the central surface brightness 
derived from this method and the bulge-disk decomposition. 

The central surface brightness is then corrected for the ge- 
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Fig. 1.— Distributions of stellar mass, gas mass, stellar-to-gas mass ratio and surface brightness of the spectroscopic sample (top) and the photometric sample 
(bottom). 


ometric effect of the inclination as follow; 


Mo,; = /^o - 2.5 log(Z7/a), (2) 

where a and b are the major and minor axis of the galaxy, 
respectively. The inclination-corrected central surface bright¬ 
ness, ^ 0 ,;. is the quantity we use for analysis in this paper. 


2.3. Stellar and Gas Masses 


For stellar masses, we adopt the color -dependent mass-to - 
light ratio from W 1 and W2 magnitudes ( |Cluver et al.|2014] l: 

logT = -1.93(Wlvega- W2vega) ”0.04, (3) 


where Wlvega = VT1ab-2.699, W2vega = W2ab-3339, and T 
is the mass-to-light ratio atWl band. T is calculated for every 
annulus until the largest annulus with reliable W2 photometry 
and then fixed for larger radius. Our sample galaxies have an 
typical mass-to-light ratio T = 0.57. 

The H I mass is calculated as Mhi = 2.356 x 10^ x x F, 
where D is distance in Mpc and F is H I flux in Jy km s“^. 
We adopt Matom = 1.4 x Mhi to include the contribution from 
helium and metals. The H I flux for the spectroscopic sam- 
ple comes from t he Extragalactic Distance Databasfl(EDD; 
Tully et al.|2009| l. Eor the photometric sample, we Mopt the 
HI flux from ALFAEEA. 

We do not have direct observations of the molecular gas 
towards our galaxies. To estimate the molecular gas con¬ 


tent of our galaxies, we refer to the result of |Bothwell et al.| 
(|2014b, who measured CO gas content towards ~ 100 galax¬ 


ies. [Bothwell et al.] ( |2014[ ) found that, the molecular-to-atomic 
gas mass ratio {Mh^Mhi) has a positive correlation wit h stel- 
lar mass of galaxies with some scatter. Erom figure 2 of |Both-| 
[well et aT] ( 2014|l, we estimate the molecular-to-atomic gas 


mass ratio: 


^ http://edd.ifa.hawaii.edu/ 


log(MH2/MHi) = 0.66 X log(M*/Mo)-7.392 (4) 

Therefore, the molecular mass is estimated from the com¬ 
bination of stellar and H I mass. Throughout this paper, ”gas 
mass” refers to the total gas mass, Matom+Mni- We wi ll dis - 
cuss the uncertainty and possible systematics in Section [5^ 
Eor distances, we use EDD as the main source. We adopt 
the group distance, Dgp, in the “Cosmicflows-2 Distances" 
section of EDD as our primary distance. For distances not 
available in EDD, we take distances from the ALEALEA sur¬ 
vey and converted to FIq = 75 km s“*Mpc“' in order to be con¬ 
sistent with EDD distances. 


3. THE SURFACE BRIGHTNESS DEPENDENCE ON METALLICITY 

In this section, we present the correlation between metal- 
licity and surface brightness from two methods in two inde¬ 
pendent samples. Going one step farther than existing studies, 
our analysis considers both the stellar mass and gas mass, and 
shows that the surface brightness dependence does not come 
from correlations between surface brightness and stellar mass 
or gas mass alone. 


3.1. Gas-phase Metallicity from Spectroscopy 

All galaxies in our spectroscopic sample ha ve central metal- 
licity and metallicity gradient reported in |Pilyugin et al. 
(j2014a]l. The first step is to verify whether the result in pre- 
vious studies, that the higher surface brightness galaxies have 
higher metallicities at given stellar mass, still holds if we take 
the average metallicity of the whole galaxy instead of the 
central metallicity as measured from fiber spectra. However, 
without knowing the spatial distribution of the gas, we cannot 
directly convert this information into average metallicity over 
the whole galaxy. Here we estimate the average metallicity by 
assuming an exponential gas disk, whose scale length scales 







































4 


with the size of the optical disk (R 25 ) (Bigiel & Blitz|20r2)l: 


-‘gas 


„(- 1 . 65 xr/fi 25 ) 


We calculate the average metallicity out to R 2 s: 

f^25 

(Z' 


/q “ Tigas (Zo + dZ/dr X r) lirrdr 


Ic 


^25 


Y.gas2T:rdr 


( 5 ) 


(6) 


where (Z), Zq, and dZ/drds& the average metallicity, the cen¬ 
tral metallicity and the metallicity gradient, respectively. 

Figure shows the central metallicity as a function of 
the stellar mass of the spectroscopic sample. The sample is 
divided into two subsample by surface brightness. First of 
all, the central metallicity shows the mass-metallicity rela¬ 
tion as expected. Overall, higher surface brightness galaxies 
are more metal-rich, mainly due to their more massive na¬ 
ture. But it is also clear that, at M* < IO'^'^Mq, most higher 
surface brightness galaxies have the highest metallicity for 
given mass, while no obvious distinction exists between two 
subsamples at higher mass. This is qualitatively in agree¬ 
ment with previous studies by SDSS fiber spectra ([Ellison] 
et al.||20()8| [Liang et al. [2OT0 Salim et ar][2014[ l, where for 
a given mass, high surface brightness, or compact, galaxies 
have higher metallicities, and the distinction is larger at the 
low mass end. 

Next, we show the average metallicity in Figure |^. Here, 
at M* < 1 O'°'^M 0 , the distinction between two samples is 
not as clear as in Figure |^, but galaxies in the higher sur¬ 
face brightness subsample still tend to locate at the top half of 
the mass-metallicity relation. We have computed the average 
metallicity with a few different assumptions of gas profiles 
(see Section for profiles used) and the results are similar. 
By com paring Figure 2^ and Figure]^, we found that, as ar¬ 
gued by [Ellison et al. p008| l, the combination of abundance 
gradient and aperture has its effect, but does not account for 
all the surface brightness dependence. Using the same data 
set, [Pilyugin et al.[(|2014b|) examined the relation between the 
disk surface brightness and the metallicity at the center and 
R 25 - They found that either at the center or R 2 S, galaxies with 
higher surface brightness on average have higher metallici¬ 
ties. This result also supports that the surface brightness of 
disk has an effect on the metallicity. 

The gas-phase metallicity is the ratio between metals and 
hydrogen. The nominator, the metal content, is related to the 
stellar mass, where the denominator is the total gas content 
of the galaxy. In observation, metallicity has been found to 
be tightly c orrelated with both the gas mass and stellar - to-gas 
ratio (Bothwell et al.[|2OT3 Hughes et al.[[2013j [Zahid 


et al.[|2014) l. If lower surface brightness galaxies are on av¬ 

erage more gas-rich, it would provide a direct explanation of 
the observed surface brightness dependence. A proper com¬ 
parison will be between galaxies with both same stellar and 
gas masses but difference surface brightnesses. 

To make such a comparison, we pick out pairs of galax¬ 
ies with similar stellar and gas mass {dM < 0.3 dex) from 
the sample and compute the difference in surface brightness 
(A/io,;) and metallicity {/^{O/H^avg) between the two galax¬ 
ies. Eigurej^shows the metallicity deficiency of a lower sur¬ 
face brightness galaxy compared to a higher surface bright¬ 
ness galaxy with similar masses. Each data point represents a 
comparison between two galaxies. A galaxy can appear more 
than once in the figure because there can be multiple galaxies 
with similar stellar and gas mass. We divide the data points 


into 10 equal-size Ap,o,,' bins and show the median and the 
16th and 84th percentile of the distribution in each bin. 

The difference, A[(9/i/], starts from 0 at A/io,, = 0, and 
gradually becomes more negative as Apo,; increases (lower 
surface brightness). This result indicates that, at the same stel¬ 
lar and gas mass, the metallicity depends on surface bright¬ 
ness of galaxies. The difference between low and high sur¬ 
face brightness galaxies seen in EigureS cannot be entirely 
contributed to their gas content. We further split the sample 
into two stellar mass bins. The surface brightness dependence 
is more prominent in the lower mass bin. The surface bright¬ 
ness effect is presented only at A/tq ,,•>2 in the higher mass 
bin, therefore it is not identified in Eigurej^ 


3.2. Stellar Metallicity from Broadband Colors 


We have seen that the average gas-phase metallicity of 
a galaxy depends on its surface brightness, and the effect 
is more prominent at lower masses. It is intriguing to see 
whether the same dependence also presents for stellar metal¬ 
licity. Without metallicity measurements from spectroscopy, 
we use broadband color as a proxy for stellar metallicity. The 
broadband color of an unresolved population is affected by 
both stellar age and metallicity. Therefore inferring metallic¬ 
ity from broadband colors alone is non-trivial. Previous stud¬ 
ies showed that, instead of optical colors alone, combining 
optic a l and NIR colors could break the age-metallicity degen¬ 
eracy ( [Galaz et al.[2002||Mac Arthur et al.[2004[[Eminian et ah] 
[2008|l . Figure | 4 | shows stellar population synthesis models of 
various metallicities and star formation histories generated by 


the stellar population synthesis code FSPS (Conroy & Gunn| 


2010 1 on the z' -W\ v.s. u'- g' color-color diagram. The 
zl -W 1 color is sensitive to metallicity and insensitive to age, 
except for galaxies dominated by young stellar populations. 
At u! - g' > 0.9, at fixed metallicity, the change in z' -W\ 
color with age is less then 0.1 magnitudes. Even extending to 
the youngest population, the z'-W\ changes by at most ~ 0.2 
mag. Therefore, the z'-WI color serves as a good proxy for 
stellar metallicity. Although different models can yield dif¬ 
ferent metalliciti es from broadband colors as pointed out in 
previous studies ( [Lee et al.|2007| [Eminian et al.[|2008p , here 
we are only interested in relative metallicity changes but not 
the actual value. The model dependence is then less impor¬ 
tant. 

Eigurej^shows z!-Wl colors of the photometric sample as 
a function of stellar mass. Similar to Eigure we split the 
sample into lower and higher surface brightness subsamples. 
Here we include only galaxies with u' -g' > 0.9 to avoid the 
effect from young populations. We plot the median z' -W\ 
colors of each subsample every 0.5 dex, for ^ ^ 

lO'^'^M©. 

Eirst, overall, more massive galaxies generally have larger 
z!-W\ color, indicating higher stellar metallicity. This con- 
clusion is in line with the result from SDSS fiber spectra ([Gal-[ 
lazzi et al. 2005[l. Moreover, similar to the gas-phase metallic¬ 


ity, higher surface brightness galaxies also have larger z!-WI 
colors at given stellar mass. In all 5 stellar mass bins, the high 
surface brightness subsample is ~ 0.15 mags redder than the 
low surface brightness subsample. Here we examin e the effect 
of surface brightness averaging the whole galaxy. [Gonzalez [ 
Delgado et al.[(2014[l discussed the effect of local surface 


densities on local stellar metallicities using lEU data. They 
concluded that the local surface brightness regulates metallic¬ 
ity, where denser regions have on average higher metallicities. 
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Fig. 2.— Left: Central metallicity as a function of stellar mass, color-coded by surface brightness. The higher and lower surface brightness galaxies are red 
triangles and blue circles, respectively. The surface brightness dependence is more prominent at low mass end. Right: Average metallicity as a function of stellar 
mass. Higher surface brightness galaxies are mostly located above the mean at low mass, but the difference is not as prominent. 
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Fig. 3.— Comparisons between galaxies with similar stellar and gas 
mass. Each data point comes from two galaxies with Alog(M*) < 0.3 and 
Alog(Mj) < 0.3. The horizontal axis is the difference in surface bright¬ 
ness, calculated from the lower surface brightness one (larger /io./) minus 
the higher one (smaller fiQj). The vertical axis is the difference in average 
metallicity. We divide the whole data set into 10 equal-size A^o./ bins. The 
median and 16th and 84th percentile of the distribution in each bin is shown 
as the green en'or bar. There is a clear trend that, after taking both stellar and 
gas mass into consideration, lower surface brightness galaxies are on average 
more metal-poor. We further split the sample into 2 stellar mass bins, and plot 
the median of each subsample for each A/io./ bins in red (M* > 10 ^^'^Mq) 
and blue (M* <= IO^^'^Mq), respectively. The surface brightness effect is 
larger at low mass than in high mass. 

Our result is consistent with this study. 

Similar to Figurej^ Figure|^compares the z' -Wl colors of 
two galaxies with similar stellar and gas masses. For galax¬ 
ies with u' -g' < 0.9, we further require A(u'- g') < 0.05 
in addition to the mass limit in order to mitigate the age ef¬ 
fect. We divide the data into 15 equal-size A^o,/ bins and 
plot the median, 16th and 84th percentile of the distribution. 



Fig. 4.— Grids of stellar population synthesis models FSPS for different 
star formation histories and metallicities. The r is the decay time scale of 
exponential SFH. r = inf indicates constant star formation rate. The tf (in 
Gyr) is the time when star formation started. The model shows that the z' — 
Wl color is sensitive to metallicity but insensitive to age for older stellar 
population. At if —g' > 0.9, the z' -W\ alone gives a good estimate of stellar 
metallicity. At u' -g' < 0.9, the z' — VT1 color still can be used if the if — g' 
color is known. 


Again, the difference in color starts from 0 at A^o,; = 0 and 
becomes more negative at larger A/ro.,, indicating that lower 
surface brightness galaxies have systematically lower stellar 
metallicities for a given stellar and gas mass. 

We have estimated the internal extinct ion following the pre¬ 
scription of|Verheijen & Sancisi|(|2001)l based on the H I line 
width and the inclination of galaxies, which we have the in¬ 
formation for all sample galaxies. However, this prescription 
does not consider the surface brightness of galaxies. If the 
metallicity of galaxies depends on the surface brightness, the 
internal extinction would depen d on the surface brightness as 
well. Thus, the prescription of |Verheijen & Sancisi| (|2001 1 
could systematically overestimate the internal extinction of 
low surface brightness galaxies if they are metal-poorer as 
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Fig. 5.— The z' — Wl color as a function of stellar mass of the photometric 
sample. Only galaxies with u' —g' > 0.9 are shown. The higher and lower 
surface brightness galaxies are red triangles and blue circles, respectively. 
We plot the median z' - Wl colors in every 0.5 dex M* bin between 8 < 
log(Mt/M q) < 10.5, where each subsample contains at least 15 galaxies. 
The z' — Wl color increases as stellar mass increases as expected. Lower 
surface brightness galaxies on average have more negative z' — Wl colors 
than higher surface brightness galaxies in all stellar mass bins shown. 



Fig. 6. — The comparison between galaxies with similar stellar and gas 
masses. Each data point comes from a pair of galaxies with Alog(Mt) < 0.3 
and Alog(M^) < 0.3. For galaxies with u' - g' < 0.9, we further require 
A(u' —g')< 0.05. The data set is divided into 15 equal-size /ro.i bins and the 
median and 16th and 84th percentile in each bin is shown by the green eiTor 
bar. At the same stellar and gas masses, lower surface brightness galaxies 
have smaller z' — Wl colors, indicating lower stellar metallicities. 


suggested in previous studies. An arguably better method is to 
estimate from the flux ratio between total infrared (TIR) and 
UV luminosities, i.e., the attenuated and un-attenuated light 


(Kong et ah] 2004^puat et al.|[2005[ |Gil de Paz et al. 


However, only ~ 20% of our galaxies are detected in the W4 

band, so this method cannot be applied in our case. When IR 
data are not available, a common alternative is using UV spec¬ 
tral slope or UV color to infer the TIR-UV flux ratio ( |Cortese| 
|et al.|20d6) [Salim et al.|2007| l because there is a positive cor- 
relation between the TlR-UV flux ratio and the slope of the 
UV spectrum (t he IRX - (3 relation; Kong et al.|2004[ Gil de 


jPaz et al.|2007|l. But meanwhiles, the correlation is found to 


depend on ages and dust properties of galaxies ([Kong et al. 
|2004l|C^se et al.|2006[|Gil de Paz et al.|2007| l. 


To make sure our re sults are not inherited from the poten¬ 
tial systematics of the jVerheijen & Sancisij ([2001 jl prescrip¬ 
tion, we also carry out the analysis with internal extinction 
estimated from the observed UV colors for galaxies with both 
GALEX FUV and NUV i maging (^ 90% o f our samples). 
We follow the calibration of jHao et al.| ( |201 ij ): 


Afuv = 3.83 X [(Ut/y-Vt/y)ofo-0.022], 


(7) 


and find that both prescriptions for internal extinction give 
qualitatively similar results. 

We would like to point out that, in both our samples, the 
surface brightness dependence is only moderate comparing 
to variations from individual galaxies. The individual varia¬ 
tion may dominate in galaxy samples which are too small in 
size or do not cover a wide enough range of surface bright¬ 
nesses. Therefore the surface brightness effect may not be 
detected. It will be possible in the near future to draw a 
cleaner picture with up coming integral fiel d unit (IF U) sur- 
veys such as MANGA ([Bundy et al.||2015|l or SAMI ([Brynt| 
et al.|2015] l, which can provide spatially resolved metallicity 
maps of large number of galaxies over wide ranges of mass 
and surface brightness. 

Also, we note that under the assumption of an exponential 
disk our result on the surface brightness dependence could be 
interpreted as a scale length (h) dependence, where galaxies 
at the same mass with higher central surface brightness have 
smaller scale lengths. 

Also, we note that under the exponential disk assumption 
our result on the surface brightness dependence could be in¬ 
terpreted as a scale length {h) dependence, because galaxies at 
the same stellar mass with higher central disk surface bright¬ 
nesses have smaller scale lengths. In general, galaxies are not 
perfect exponential disks, and this conversion between po and 
h will depend on the structure of galaxies. For our cases, we 
have repeated the similar analysis shown in Figure|^and Fig¬ 
ure]^ but as a function of Ah, and find similar results. 

4. THEORETICAL CONSIDERATIONS 

A large number of studies have modeled the observed mass- 
metallicity relation with different approaches and formula¬ 
tions. It is widely agreed that the close-box assumption yields 
too high metallicity. Outflows removing metal-enriched ISM 
and/or inflows of pristine gas diluting ISM are required to 
quantitatively match the observed mass-metall i city relation 
([Spitoni et al.[[2010[ [Peeples & Shankar[[201 1| [Zahid et al.| 
[2012[ [Lilly et al.|2013[[Zahid et al |2014| l. 

While in their work, the spatially integrated mean relation 
is successfully reproduced, there is little work addressing the 
scatter around the mean. We have shown in Section that, 
for a given mass, the metallicity is correlated with the surface 
brightness. It is intriguing to investigate how the spatially re¬ 
solved surface brightness and metallicity profiles are affected 
by the inflow and outflow properties. 

For this purpose, we apply the chemical evolution model of 
Kudritzki et al.[ ([2015J, which addresses the relation among 
inflow, outflow, metallicity and stellar-to-mass ratio in the 
spatially resolved case. The inflow and outflow of a galaxy 
are parameterized by two factors: 


and 


^ _ Alf^ccr 

-tP 

A^loss 


v = 


V' 


( 8 ) 


(9) 

































































Surface Brightness and Mass-metallicity Relation 


7 


where Mioss and Maccr are the mass-loss and mass-accretion 
rate, and ijj is the SFR. 

This model makes a few assumptions. First, ry and A are 
assumed to be constant in time. This assumption constrains 
only the ratio between the SFR and inflow and outflow rate, 
therefore, all three quantities can be time variables. Second, 
the inflow gas is free of metals, and the outflow gas has the 
same metallicity as the ISM at the time launched. Third, the 
nucleosynthetic yield (y) and the stellar mass return fraction 
(R) are both constant. Under these assumptions, the metallic¬ 
ity can be analytically expressed as 


Z(t): 


1 - 


i-r’mm) 




( 10 ) 


Therefore, the metallicity is determined by p. A, and the 
stellar-to-gas mass fraction at a given time. If a galaxy has 
spatially resolved information on stellar mass, gas mass, and 
metallicity, the three constraints can be used to And the best- 


fit rj an d A of the galaxy. We refer readers to Kudritzki et al. 
( |2015[ ) for detailed derivation and discussion of the model. 

We apply the model only to the spectroscopic sample but 
not the photometric sample. The stellar metallicity inferred 
from broadband colors can depend on stellar populati on syn¬ 
thesis models and is less certa in for quantitative work ( |Lee et| 
|al.|2007]|Eminian et al.|2008| . 

For each galaxy in the spectroscopic sample, we have both 
the metallicity and stellar mass profiles from observations. 
For the gas distribution, we follow the form in Equation!^ and 
assume an pure exponential gas disk with slope scaled^with 
R 25 - This assumption is supported by observations between 
0.2 and 1.0 R 2 S, but is not valid for radii outside this range 
(|Bigiel & Blitz|2012|). Therefore, our nominal range for anal- 
ysis is limited between 0.2 R 25 and 1.0 ^ 25 - For the amount 
of gas within this range, we do the analysis under two sep¬ 
arate assumptions. The first one is distributing the observed 
amount of gas exponentially from galaxy center to infinity. 
Therefore, the gas surface density is determined by both R 25 
and the integrated observed gas content. But the validity of 
this assumption is potentially affected by the gas profile be¬ 
yond R 25 - Thus, we assume as an alternative a g as surface 
density profile, which is only constrained by R 25 dBigiel & 
|Blitz|20121 l: 


= 2.1 X exp(-1.65 X r/R 25 ),'Stmn ^ 14M©/pc^. (11) 


Although the second assumption, determining gas density 
profile purely from the size of stellar disk, may see m to be 
less convincing, it is in fact supported by observations (|Bigiel 
& Blitz||2012| l. Comparing the two assumptions with the 19 
galaxies studied in [Kudritzki et al.| ( |2015| ), which have a re¬ 
solved gas distribution, we find that the first assumption on 
average yields ^0.17 dex lower gas surface density, while 
the second assumption under-predicts gas surface density by 
^ 0.11 dex within the nominal radius range. As we will And 
out later, this range of uncertainty is not a concern for our 
discussion. 

For each galaxy we have now the photometrically derived 
stellar mass column density profiles and gas column density 
profiles either described by Equation|^and the integrated gas 
mass or by Equation!^ We now calculate model metallicity 
profiles from the stellar^to-gas mass ratio profile and pairs of 
77 and A between 0 and 3 for every hundredth interval. We 


adopt R = 0.4 and y = 0.00313 as calibrated in Kudritzki etj 
ari(|2015|l. A 0.15 dex is added to the H Il-region metallicity 
measured by Pilyugin et al. ( 2014a 1 to match the metallicity 
zero point obtained by [Kudritzki et al.] ( [2015| l as the result of 
stellar spectroscopy in spiral galaxies. We then compare the 
model metallicity profile with the observed metallicity profile 
and compute the from the difference between two profiles 
and the uncertainties for every 0.05 R 25 from 0.2 R 25 to 1.0 
R 25 or the last reliable measurement of surface brightness for 
adopted pairs of rj and A. The pair of p and A with the mini¬ 
mum is then defined as the best fit model. EigureH shows 
an example of our comparison between the model andthe ob¬ 
servation. We note that the assumption on gas profile is valid 
only in a statistical sense. There is ^ 60% scatter among in¬ 
dividual galaxies. Therefore, we will only discuss the average 
properties from the model. 



Fig. 7.— An example of the model fit. The black line represents the ob¬ 
served radial metallicity profile, with the la uncertainty indicated by the 
dashed lines. The best-fit model metallicity profile is plotted as the blue line. 
The close-box model (77 = 0, A = 0) is shown by the red dash-dotted line as a 
reference. 


Our result with two different gas profiles is presented in 
Eigure|^ Galaxies are first sorted into 3 bins in baryonic (Eig- 
ure|^ and c) or stellar mass (Eigurej^ and d). In each mass 
bin, galaxies are further separated into two equal-size bins by 
their surface brightness. Each bin contains ^ 20 galaxies. We 
show the median and the 16 and 84 percentile of the distri¬ 
butions of rj and A in each sub-sample. The difference be¬ 
tween the median surface brightnesses of two sub-sample is 
^ 1 mag for the high and mid-mass bins, and ~ 2 mags for 
the low-mass bin. As a comparison, we also plot the best-fit rj 
and A from the 19 galaxies with measure d stellar-to-gas mass 
profiles studied by [Kudritzki et al.[ ( [2015[ ). 

To investigate the uncertainty resulting from the parameter¬ 
ized gas profile, we considered four more different gas pro¬ 
files. Eollowing the la scatter of the slop e of the universal 
gas profile given in [Bigiel & Blitz[ ([20T2|, we first consider 
two cases where the gradient of the gas distribution is -1.78 
or -1.49, instead of the canonical value of -1.65. Eor the other 
two cases, we Ax the gradient, but double or half the amount 
of gas in order to access the uncertainty from the assumed 
amount of g as and the effect from possible systematics (see 
Section [572| for more discussion). Eor all four cases, the me¬ 
dian p and A of each sub-sample is calculated and plotted in 
lower panels of Eigurej^ Variation in slope is shown by light 
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Fig. 8.— Best-fit rj and A as a function of mass and surface brightness, (a): gas amount scaled by R 25 , binned by baryonic mass, (b): gas amount scaled by 
R 25 , binned by stellar mass, (c): gas amount scaled by S 25 measured gas mass, binned by baryonic mass, (d): gas amount scaled by R 25 and measured gas 
mass, binned by stellar mass. Galaxies in the spectroscopic sample are first sorted into 3 mass bins then two surface brightness bins in each mass bin. Higher 
surface brightness galaxies are in red and lower surface brightness galaxies are in blue. In the upper panel of each subplot, the error bars represent for the 16th 
and 84th percentile of the distribution of r} and A in each subsample. The filled symbols in lower panel of each sub^t show the uncertainties fro m varying the 
slope of surface density of gas. Open symbols illustrate the uncertainties from varying the gas amount (see Section|4). The 19 galaxies studied by|Kudritzld et| 
|al.|j2015} are also plotted as a comparison. 


filled symbols and variation in gas amount is shown by open 
symbols. 

Generally, both p and A depend on the mass of galaxy, 
where higher mass galaxies have on average lower 77 and A. 
However, rj depends strongly on the gas amount assumed. 
Without knowing the amount of gas, p is essentially not con¬ 
strained in the lowest mass and lower surface brightness bin. 
In all subsamples, the difference in rj between high and low 
surface brightness bins could come from systematically dif¬ 
ferent assumed gas amount, if there is any. On the other hand, 
A is more robust against variations in gas profiles. In the two 


low mass bins, lower surface brightness galaxies have larger 
A regardless of gas profile assumed. 

The same result also can be seen in galaxies with observed 
stellar-to-gas mass profiles (open diamonds). Only A in the 
low mass bin shows a surface brightness dependence, while 
rj, and A in the high mass bin does not seem to depend on 
surface brightness. However we note that this result comes 
from a small sample of only 5 galaxies in each bin. 


5. DISCUSSION 
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5.1. Implications from the Mass and Surface Brightness 
Dependence on A 

Chemical evolution model suggests that, for given mass, 
lower surface brightness galaxies tend to have larger A. In this 
section we discuss possible explanations for the correlation 
between A and the surface brightness of a galaxy. 


5.1.1. Star Formation Efficiency 

The accretion parameter A is the gas accretion rate divided 
by star-formation rate. It can be understood as the recipro¬ 
cal of star formation efficiency; while gas is falling onto the 
galaxy, A monitors how much of it is turned into stars. A 
larger A indicates inefficient star formation. 

Our model suggests that lower surface brightness galax¬ 
ies are less efficient in turning gas into stars, which is in 
agreement with the prevalent idea. Observations of giant 
low surface brightness galaxies (M* > IO'^Mq) showed that 
they have copious amou nts of H I gas that is comparable to 
their stellar components ( [Matthews et al.|2001[ [O’Neil et al.j 
2004|l. Their gas-to-stellar mass ratio is higher than their 
high surface brightness counterparts ( Burkholder et al.|2001[ 


O’Neil et al.|20()7|l. By com paring star formation rate and HI 
mass,[Boissier et ( 2007[) found that low surface brightness 


galaxies have lower SFR compared to high surface bright¬ 
ness galaxies with similar HI mass, evidence that low surface 
brightness galaxies are inefficient in star formation. 

The low star formation efficiency in low surface brightness 
galaxies is also suggested by theory that low surface bright¬ 
ness galaxies are stable against both global and local growth 
of instability ([Dalcanton et al.[1997[[Mihos et al.[1997[ Mayer[ 


[& Wadsley et al.|2004[|Cjhosh & Jog[2014[l. Stability analysis 

on individual low surface brightness gal axies demonstrated 
that they have high Toomre Q parameter ( Mihos et al.|1997 
Ghosh & Jogj2014 1 , where Q = a k/3.36 G II, > 1 for stabil¬ 
ity, and a and k are the radial velocity di spersion of the stars 
and the epicyclic frequency, respectively ( [Toomre |1964[ l. 

Structures like spiral arms and ba rs are also harder to form 
due to a lack of global instability ( [Mihos et al.[[1997 l. Gas 
clouds in low surface brightness disks experience less cloud 
collisions and compression, which can trigger star forma¬ 
tion. This proposition is also supported by the fa ct that fewer 
bars a re seen in low surface brightness galaxies ( [Impey et al.[ 
1996|l. The larger A inferred from our model, thus, agrees 


with the majority of previous studies in terms of star forma¬ 
tion efficiency. 


5 . 1 . 2 . Environmental Effect 

Since A is the ratio between mass accretion rate and star 
formation rate, for two galaxies with the same stellar mass, 
the one with higher A, i.e. lower surface brightness, should 
have had more inflow material over time. This may imply that 
low surface brightness galaxies reside in higher density region 
than higher surface brightness galaxies of the same mass, or 
more likely to be the central galaxy in a halo because central 
galaxies have overall hi gher accretion rat es shown by hydro- 
dynamical simulations ( [Keres et al.|2009| l. 

This assertion seems to be at odds with the general percep¬ 
tion, that low surface brightness galaxies have fewe r neigh¬ 
boring galaxies on scales of < 1 Mpc to several Mpc (Bothun 


[et al.[ 1993[ [Rosenbaum et aL|[20()9t [Galaz et al.[[2011[ l. But 

this is not necessary in conflict with our results. While there 
are several ways to define the surface brightness of a galaxy, 
low surface brightness galaxies in the literature refer mostly 


to galaxies whose surface brightness are fainter than a cer¬ 
tain value. This definition yields a heterogeneous sample. A 
galaxy with low surface brightness could be a result of an ex¬ 
tended distribution of baryons, its less massive nature, or be¬ 
ing massive but gas-rich (star-poor) condition. Each type of 
galaxies may have its own environmental dependence, and the 
c ontribution from each typ e is not clear. 

Rosenbaum et al.[(2009|l found that low surface brightness 


galaxies reside in underdensed regions at scales > 2 Mpc from 
SDSS DR 4 sample. In the same study, they further split 
the sample into two redshift bins and showed that this en¬ 
vironmental dependence is more p rominent in the low-z bin 
than in the high-z bin. As shown in[Rosenbaum et al.[([2009[), 
their low-z low surface brightness sample consists mainly of 
dwarf galaxies, and the high-z low surface brightness sam¬ 
ple is dominated by regular spirals. The behavior in two red- 
shift bins may suggest that, at least part of the apparent en¬ 
vironmental effect is inherited from the luminosity, or mass, 
of the galaxy. In Figure we are comparing galaxies with 
the same mass, where the difference should re flect the effect 
from g alaxy structure but less from mass. Also, [Blanton et ah 


( [2005 1 studied the relationship between environment and op¬ 
tical properties of galaxies in the SDSS and found that, lower 
surface brightness galaxies are in denser region for a given 
luminosity. The result from[Blanton et al.[([2005|l and our re¬ 
sult suggest that it would be necessary to investigate each type 
of galaxies separately and not put heterogeneous individuals 
under the same umbrella. 

Although our result indirectly suggests that galaxies with 
more extended disks tend to reside in denser regions, we real¬ 
ize at the same time that the question whether the environment 
affects galaxy structure is beyond the scope of this paper and 
worth its own dedicated study. 


5 . 1 . 3 . Assembly History 

The discussion so far is under the assumption of constant 
A. This assumption has the physical motivation that inflow 
increases a gas reservoir therefore the capability of forming 
stars, but it is hard to verify observationally. Some observa¬ 
tions o f high-redshift massive galaxies (1 < z < 3) suggested 
A > 1 ( jTacconi et al.|201^[2013[|Yabe et al.|2015| l, which is 
larger than what we derive for local galaxies (A < 1) there¬ 
fore constitutes a challenge to our assumption of constant A. 
Nevertheless, it is not clear whether galaxies in different stud¬ 
ies can be directly compared, i.e., our sample is not necessary 
the descendant of galaxies in previous studies. Moreover, in 
spite of possible flaws in the assumption, this model has been 
shown to be able to reprod uce the distribut ion of metallicity 
gradients of local galaxies ( [Ho et al^[2015[ l. Practically, the 
metallicity is not sensitive to inflow and outflow in the gas- 
rich regime. Because metallicity is the content of metals nor¬ 
malized by the amount of gas, inflow and outflow would not 
drastically change it when the galaxy still possess a vast gas 
reservoir. Therefore, the model parameters, rj and A, reflect 
the more recent accretion activities and less the early phase of 
galaxy formation. Although the assumption of constant rj and 
A may not necessarily hold, we consider our fitting parame¬ 
ters should at least reflect the inflow and outflow properties at 
recent times. 

If we assume that A is decreasing over time, as hinted 
by observations, the dependences on surface brightness we 
see from local galaxies are qualitatively consistent with the 
canonical view of halo formation. More compact galaxies are 
considered to reside in compact halos with lower spin which 
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form relative early ( [Jimenez et al.|1998 Mo et al.|1998 1 . For 
baryons, hydrodynamic simulation showed that most baryons 
in low spin halos co llapse into the cen ter of the galaxies within 
the first ~ 1.5 Gyrs ( [Kim & Lee|2013| ), therefore implying low 
accretion at the present day. 

Overall, galaxies with higher surface brightness have likely 
passed their peak accretion phase, while low mass, lower sur¬ 
face brightness galaxies are still assembling their baryons. 
The low A observed in the high mass bin could also re¬ 
sult from the virial shocks developing in halos with Mhaio ^ 
IO' ^Mq, which suppress co l d material falling o nto the galax¬ 
ies ( |K^es et al.|2005||2009| [Brooks et al.|2009| ). In this case, 
the mass is the dominant factor, and the surface brightness ef¬ 
fect is relatively small and, therefore, not observed. Both fac¬ 
tors could lower the accretion rate, Maccr in high mass, high 
surface brightness galaxies. On the other hand, it has been 
shown that low surface brightness galaxies on av erage have 
lower SFR ([O’Neil et al.|2007[[Huang et al.|2012|l, which af¬ 
fects the denominator of A. ITierefore, low mass, low surface 
brightness galaxies would have an average higher A, from the 
combination of stronger accretion and lower SFR. 


5.2. Possible Systematics 
5.2.1. Molecular gas content 


In Section!^ we compare galaxies with the same stellar and 
gas mass, where molecular mass is estimated from a stellar- 
mass-dependent M(H 2 )/M(HI) fraction from Bothwell et al. 
([2014[|. Under this assumption, galaxies with the same stel¬ 
lar and HI mass will have the same amount of molecular gas. 
However, one could naturally expect that lower surface bright¬ 
ness galaxies have a deficit in molecular gas due to their lower 
metallicity, so that the neutral gas will not be able to cool 
enough to form H 2 . 

Detections of molecular gas in low surface brightness 
galaxies suggested that they have low M(H 2 )/M(HI) fraction 
compared to galaxies with higher surface brightness (fO’Neil 


jet al.|2003||2004l|Matthews et al.|2005l|Das et al.|20(^6|l. As 


a result, our estimate could lead to systematic overestimate of 
total gas content in lowe r surface brightness ga laxies. 

Based on the result of [Bothwell et al.| ( |2014| ), the scatter of 
the M{H 2 )/M{HI) fraction around the mean value is ~ 0.4 dex 
for given stellar mass. Some of the scatter is possibly inher¬ 
ited from the correlation with metallicity and surface bright¬ 
ness. According to Equation]^ the H 2 mass is less than 10% 
of H I at M, < 1 O'°M 0 . The 0.4 dex scatter corresponds to 
< 25%, or < 0.1 dex, of total gas mass, and is therefore neg¬ 
ligible. 

Of more concern is the high mass end, where M* ~ 
1 O"M 0 . At this stellar mass, the average H 2 mass is compara¬ 
ble to HI mass, and the scatter in total gas mass is ~ 0.25 dex. 
To account for this possible systematic effect, we have al¬ 
lowed a 0.3 dex difference in mass for galaxies to be con¬ 
sidered as ’’the same mass” in Section Moreover, for our 
chemical evolution model, we have run the test cases of differ¬ 
ent gas content in Section|^and shown the result in Figure]^ 
Our discussion on A is not affected by changing the gas con¬ 
tent within the nominal range. Especially at low mass, where 
the surface brightness effect is more prominent, the molecular 
gas content is essentially negligible. 


5.2.2. Metallicity scale 

In Section]^ we a dd a correction of 0.15 dex to the metal- 
licities measured by|Pilyugin et al.|([2014a|l. This correction 


has been found by Kudritzki et al. ( 2015(1 and is the result of 
a comparison of nietallicities obtained from detailed quanti¬ 
tative spectroscopy of blu e supergiant stars with H Il-region 
metallicities measured by |Pilyugin et al.| ( [2014a[ l. The yield 
(y = 0.00313) used in our chemical evolutio n model of sec- 
tion 4 has been empirically determined by [Kudritzki et al.[ 
([2015[l from the metallicity and metallicity gradient in the 
IVlilky Way disk as observed from Cepheids and B-stars The 
yield and the metallicity scale used in Section|2are,thus, con¬ 
sistently based on quantitative stellar spectroscopy. 

The inflow and outflow strength determined in our model 
could be affected by systematic uncertainties of the metallic¬ 
ity scale use d. To access this uncert ainty we apply a similar 
procedure as [Kudritzki et al.| ([2015[l We repeat the analysis 
of Section [4 but this time adding corrections to the Pilyu-[ 
[gin et al.] ( [2014a[ ) metallicities of 0.10, 0.05, and 0.00 dex, 
respectively. Figure shows the effects from changes of the 
zero point. While the zero point is important in a quantitative 
sense (decreasing the metallicity leads to larger values of 77 
and A), it does not affect our major conclusion. Lower sur¬ 
face brightness galaxies have higher A values than higher sur¬ 
face brightness galaxies of similar mass. In Eigure|^we only 
show the results for one type of gas profiles and only binned 
by baryonic mass corresponding to the case in Eigure|^. Ap¬ 
plication of the zero point shifts to the other cases lead to a 
similar conclusion. 

We note that changing the zero-point of the metallicity scale 
is equivalent to changing the yield in the chemical evolution 
model, as can be seen in Equation [T0| The result discussed 
here can also be seen as a test of the mnuence from the uncer¬ 
tainty of the yield. 


6. SUMMARY 

In this paper, we investigate whether the metallicity of 
galaxies depends on their surface brightnesses and discuss the 
possible origin of the dependence. 

1. Previous studies found that, at given stellar mass, low 
surface brightness galaxies are more metal-poor. These 
studies used SDSS fiber spectra, where the aperture ef¬ 
fect is difficult to account for. We show that the result 
remains valid when the average metallicity of the whole 
galaxy is considered. 

2. The surface brightness dependence still exists if gas 
mass is also taken into account. Eor galaxies with simi¬ 
lar stellar and gas masses, lower surface brightness sys¬ 
tems tend to have lower metallicity. 

3. We use chemical evolution models to constrain the in¬ 
flow and outflow properties of galaxies. The ratio be¬ 
tween accretion rate and star formation rate (A) is larger 
in low mass, lower surface brightness galaxies. The sur¬ 
face brightness effect is more prominent in low mass 
galaxies. On the other hand, the outflow property is not 
well constrained. 

4. The high A in lower surface brightness galaxies can 
be directly interpreted as low star formation efficiency. 
This interpretation is consistent with the prevalent idea. 

5. Eor galaxies with the same stellar mass, lower surface 
brightness galaxies should have more inflow material 
overtime because of higher A, which is the ratio be¬ 
tween mass accretion rate and star formation rate. This 
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Fig. 9.— Test of the systematical effect from the zero point of met allicity scale on rj and A. We adopt 4 different metallicity scales, adding 0.15 dex (the scale 
adopted for this work), 0.10, 0.05, and 0.00 dex to the metallicities of|Pilyugin et al.|12014al. Gas profiles are scaled by R 25 as in Figure|8^. The zero point affects 
r] and A in a quantitative sense, but not our conclusion, where lower surface brightness galaxies have higher A values than higher surface brightness galaxies of 
similar masses. 


result indirectly suggests that lower surface brightness 
galaxies reside in denser regions than higher surface 
brightness galaxies of the same mass, or more likely to 
be the cetral galaxy in a halo, therefore they have more 
inflow material. 


6 . If A is not constant in time as assumed in the model, the 
surface brightness dependence on A may be an indica¬ 
tion of different accretion histories, where low mass, 
lower surface brightness galaxies are currently at their 
major accretion phase, while this phase in high mass, 
higher surface brightness galaxies has passed. 


In the near future, large IFU surveys will provide maps of 
metallicity, and perhaps measurements of inflow and outflow 
activities, in a large number of galaxies spanning wide ranges 
of galaxy properties. We will understand better how galaxy 
structure affects the evolution of galaxies with the aid of IFU 


data. 
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